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Abstract. We consider N single server infinite buffer queues with service rate /?. Customers arrive at rate Na, 
choose L queues uniformly, and join the shortest. We study the processes t e i-^ = (i?^(fc))fcgN for 
large N, where (k) is the fraction of queues of length at least k at time t. Laws of large numbers (LLNs) 
are known, see Vvedenskaya et al. |15|, Mitzenmacher L12J and Graham |5|. We consider certain Hilbert 
spaces with the weak topology. First, we prove a functional central limit theorem (CLT) under the a priori 
assumption that the initial data R^ satisfy the corresponding CLT. We use a compactness-uniqueness method, 
and the limit is characterized as an Ornstein-Uhlenbeck (OU) process. Then, we study the i?^ in equilibrium 
under the stability condition a < fi, and prove a functional CLT with limit the OU process in equilibrium. We 
use ergodicity and justify the inversion of limits limAr^oo limj^oo = limj^oo liniAr^oo by a compactness- 
uniqueness method. We deduce a posteriori the CLT for Rq under the invariant laws, an interesting result in 
its own right. The main tool for proving tightness of the implicitly defined invariant laws in the CLT scaling 
and ergodicity of the limit OU process is a global exponential stability result for the nonlinear dynamical 
system obtained in the functional LLN limit. 
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1 Introduction 
1.1 Preliminaries 

We consider a Markovian network constituted of > L > 1 infinite buffer single server queues. 
Customers arrive at rate Na, are each allocated L distinct queues uniformly at random, and join the 
shortest, ties being resolved uniformly. Servers work at rate f3. Arrivals, allocations, and services are 
independent. For L = 1 we have i.i.d. Mq/M^/I/oo queues. For L > 2 the interaction structure 
depends only on sampling from the empirical measure of L-tuples of queue states: in statistical 
mechanics terminology, the system is in L-body mean-field interaction. We continue the large N 
study introduced by Vvedenskaya et al. fl31 and Mitzenmacher and continued in Graham l5l . 

The process {X^)i<i<]\f is Markov, where (t) denotes the length of queue i at time t in R+. 
Its empirical measure = jf YliLi ^x^ ^'^^ samples in 7'(B(R+,N)), and its marginal process 
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= {X^)t>Q with = = jfY^Zi^x'^it) has sample paths in D(M+,T'(N)). We are 
interested in the tails of the marginals and consider 

1=1 

and {k) is the fraction of queues of length at least k at time t. For the uniform topology on 

V = {{v{k))k^„ ■ v{Q) = 1, v{k) > v{k + 1), liniT; = 0} C CO , = V n ^N^ , 

coinciding here with the product topology, the process has sample paths in B(M4., V^). 

The processes and R^ are in relation through p G 'P(N) < — > ti G V for f (A;) = p[k, oo) and 
p{A;} = v{k) — v{k + 1) for k in N. This classical homeomorphism maps the subspace of probability 
measures with finite first moment onto V n ^i, corresponding to a finite number of customers in the 
network. The symmetry structure implies that these processes are Markov. 

The stationary regime has great practical relevance. The stability condition a < [5 (Theorem 5 (a) 
in fl31 . Lemma 3.1 in fT2l . Theorem 4.2 in L5 |) is obtained from ergodicity criteria yielding little 
information. We study the large asymptotics of R^ , first for transient regimes with appropriately 
converging initial data, and then in equilibrium using an indirect approach involving ergodicity in 
well-chosen transient regimes and an inversion of limits for large N and large times. Law of large 
numbers (LLN) results are already known, and we obtain functional central limit theorems (CLTs). 

1.2 Previous results: laws of large numbers 

We relate results found in essence in Vvedenskaya et al. fTSl. We follow Graham ||5l which extends 
these results, notably by considering the empirical measures on path space jjL^ and thus yielding 
chaoticity results (asymptotic independence of queues). Chapter 3 in Mitzenmacher fT2l gives re- 
lated results. (The rates a and {5 correspond to A and 1 in ll5llT2l and v and A in 13.) 

Consider the mappings with values in Cq given for v in cq by 

F+{v){k) = a{v{k - l)^ - v{k)^) , F^{v){k) = (3{v{k) - v{k + 1)) , k>l, (1.1) 

and F = — F^, and the nonlinear differential equation ii = F{u) on V, given for t > by 

ut{k) = a{ut{k - 1)^ - ut{k)^) - P{ut{k) - ut{k + 1)) 

= aut{k - 1)^ - {aut{k)^ + I3ut{k)) - Put{k + 1) , k>l. (1.2) 

This corresponds to the systems (1.6) in lITSl . (3.5) in (12\ and (3.9) in ||5l. Note that F_ is linear. 

Theorem 1.1 There exists a unique solution u = {ut)t>o taking values in V for l\1.2i . and u is in 
C(M+, V). Ifuo is in V Cili then u takes values in V f] l\. 
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Proof. We use Theorem 3.3 and Proposition 2.3 in 13. These exploit the homeomorphism be- 
tween 'P(N) with the weak topology and V with the product topology. Then (ll.2t corresponds 
to a non-linear forward Kolmogorov equation for a pure jump process with uniformly bounded 
(time-dependent) jump rates. Uniqueness within the class of bounded measures and existence of 
a probability-measure valued solution are obtained using the total variation norm. Theorem 1 (a) in 
1151 yields existence (and uniqueness) in V H ^i. □ 

Firstly, a functional LLN for initial conditions satisfying a LLN is part of Theorem 3.4 in fSl and 
can be deduced from Theorem 2 in flSll . 

Theorem 1.2 Assume that {Rq )n>l converges in law to uq in V. Then {R^)n>l converges in law 
in I])(R+, V) to the unique solution u = {ut)t>o starting at u^for M.2\ . 

Secondly, for a < /? the limit equation (II. 2t has a globally attractive stable point u in V n ^i. 

Theorem 1.3 Let p = a/ (3 < 1. The equation M.2\ has a unique stable point in V given by 

u = {mum , m = p(^'-mL-i) = ^L-^+L-^+...+i 

and the solution u of M.2\ starting at any uq in V Ci ii is such that lim^^oo Ut = u. 

Proof. Theorem 1 (b) in (T5i yields that u is globally asymptotically stable in V n £i. A stable point 

n in V satisfies (3u{k + 1) — au{k)^ = f3u{k) — au{k — 1)^ = • • • = Pu{l) — a and converges to 
0, hence n(l) = a/f3 and n(2), n(3), ... are successively determined uniquely. □ 

Lastly, a compactness-uniqueness argument justifies the inversion of limits liiii7v_>oo limt— >oo — 
lim^^oo limiv^oo> which yields a result in equilibrium. This method, used by Whitt fl6l for the star- 
shaped loss network, is detailed in Graham j6| Sections 9.5 and 9.7.3. The following functional LLN 
in equilibrium (Theorem 4.4 in ||5l) can be deduced from (TSi but is not stated there, and impUes that 
under the invariant laws liniAr^^oo (k)) = u{k) for A; G N (Theorem 5 (c) in lITSi ). 

Theorem 1.4 Let p = a/ j3 < 1 and the networks of size N be in equilibrium. Then {R^)n>l 
converges in probability in B(M_|_, V) to u. 

Note that u{k) decays hyper-exponentially in A; for L > 2 instead of the exponential decay 
corresponding to i.i.d. queues in equilibrium (L = 1). For finite networks in equilibrium there is at 

most exponential decay since P(Xf H h > Nk) < P(Xf > k) -{ hP(X^ > A;) and 

by comparison with an Matq/Mtv/?/! queue 

E(i?f (fc)) = P(xf (t) > A;) > Ip^'^ , k>0. (1.3) 
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The asymptotic queue sizes are dramatically decreased by this simple load balancing (or resource 
pooling) procedure, which carries little overhead even for large N since L is fixed (for instance 
L = 2). This feature is quite robust and true for many systems, as was illustrated on several examples 
by Mitzenmacher IIT2I and Turner llT4l using proofs as well as simulations. It can be used as a 
guideline for designing practical networks. In contrast, the bound (II. 3t assumes the best utilization 
of the N servers, fully collaborating even for a single customer. 

Theorem 3.5 in Graham 13 gives convergence bounds on bounded time intervals [0, T] for i.i.d. 
[X^ {0))i<i<]\f using results in Graham and Meleard 0. This can be extended if the initial laws 
satisfy a priori controls, but it is not so in equilibrium (the bounds are exponentially large in T). 

1.3 The outline of this paper 

The study of the fluctuations around the functional LLN will yield for instance asymptotically tight 
confidence intervals for the process t ^ A^~^Card{z = 1, . . . , : (t) G ^4}. In a realistic set- 
ting (finite number of finite buffer queues) such confidence intervals would allow network evaluation 
or dimensioning in function of quality of service requirements on delays and overflows. The LLN 
on path space concerns objects such as A^^^Cardji = 1, . . . , : 1— > Xj^{t)) G with a richer 
temporal structure, but topological difficulties usually block the corresponding fluctuation study. 

We consider the process and solution u for (II. 2t starting at Rq in and in V, and 

= Vn{R^ -u). (1.4) 

The processes Z^ = {Z^)t>o will be studied in the Skorokhod spaces on appropriate Hilbert spaces 
with the weak topology. These spaces are not metrizable and require appropriate tightness criteria. 

We first consider a wide class of Rq and uq under the assumption that {Zq)n>i converges 
in law (for instance satisfies a CLT). We obtain a functional CLT in relation to Theorem 11.21 with 
limit given by an Ornstein-Uhlenbeck (OU) process starting at the limit of the {Zq)n>l- This 
covers without constraints on a and f3 many transient regimes with explicit initial conditions, such 
as initially empty networks, or more generally i.i.d. initial queue sizes. 

We then focus on the stationary regime for a < [3. The initial data is now implicit: the law of 
Rq is the invariant law for R^ and uq = u. We prove tightness for {Zq)n>l using the ergodicity 
of Z^ for fixed N and intricate fine studies of the long-time behavior of the nonlinear dynamics 
appearing at the large limit. The main result in this paper is a functional CLT in equilibrium for 
{Z^) ]\j>L with limit the OU process in equilibrium. This implies a CLT under the invariant laws for 
{Z^)n>l, an important result which seems difficult to obtain directly. 

Section 121 introduces without proof the main notions and results. Section|3l gives the proof of the 
functional CLT for converging initial data by compactness-uniqueness and martingale techniques. 
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We then consider uq = u. We study the OU process in Section|4l derive a spectral representation 
for the linear operator in the drift, and prove the existence of a spectral gap. A main difficulty is that 
the scalar product for which the operator is self-adjoint is too strong for the hmit dynamical system 
and the invariant laws for finite N. We consider appropriate Hilbert spaces in which the operator is 
not self-adjoint and prove exponential stability. 

In Section|5]we likewise prove that u is globally exponentially stable for the non-linear dynamical 
system. In Section |6l we obtain bounds for uniform for i > and large N, using the preceding 
stability result in order to iterate the bounds on intervals of length T. Bounds on the invariant laws 
of follow using ergodicity. The proof for the functional CLT in equilibrium follows from a 
compactness-uniqueness argument involving the functional CLT for converging initial data. 



2 The functional central limit theorems 
2.1 Preliminaries 

The exponential of a bounded linear operator is given by the usual series expansion. Let Cg and £p for 
p > 1 be the subspaces of sequences vanishing at of the classical sequence spaces cq (with limit 0) 
and £p (with summable p-th power). In matrix notation we use the canonical basis, hence sequences 
vanishing at are identified with infinite column vectors indexed by {1, 2, • • •}. The diagonal matrix 
with terms given by the sequence a is denoted by diag(a). Sequence inequalities, etc., should be 
interpreted termwise. Empty sums are equal to and empty products to 1. Constants such as K may 
vary from line to line. Let = {9^)k>i be the geometric sequence with parameter 6. 
For a sequence w = {w{k))k^i such that w{k) > we define the Hilbert spaces 



L2{w) = {x£R^ : x{0) =0, ||x| 



L2{w) 



and in matrix notation [x, y)L^{^yj) = x* d\a.g{w~^)y . We use the notation L2{w) since its elements 
will often be considered as measures identified with their densities with respect to the reference 
measure w. In this perspective Li{w) = £^ and if w is summable then ||x||i < ||^t'||i^^||2;||j;^2(u)) 
and L2{w) C £\. Using -^2(1) = ^2 ^s a pivot space, for bounded w we have the Gelfand triplet 
L2{w) (Z lie L2{wY = L2{w-^). 

Another natural perspective on L2{w) is that it is an I2 space with weights, and we consider the 
£1 space with same weights (the notation being chosen for consistency) 

= J G : x(0) = , \\x\\i^(^^) = ^ \x{k)\w{k)'^ < 00 I 

I k>l ) 

and X G L2{w) <^ x"^ e £i{w) with = The inclusion V n £i{w) ^ V n L2{w) 

is continuous since x"^ < \x\ for \x\ < 1. The following result is trivial. 
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Lemma 2.1 Ifw = 0{v) and v = 0{w) then the L2{v) and L2{'w) norms are equivalent, and the 
^l{v) and ii{w) norms are equivalent. 



In the sequel we often assume that w = {'Wk)k>i satisfies the condition that 



3c,d>0,VA;>l,0< cw{k + 1) < ^(A;) < dw{k + 1) 



(2.1) 



which is satisfied by go = {6^)k>i with c = d = 1/6* for 6* > 0. It implies that w{l)d{l/dY < 
w{k) < 'w{l)c{l/ c)^ which bounds w by geometric sequences. The norms have exponentially strong 
weights for c > 1. We give a refined existence result for ()1.2|) . (Proofs are left for later.) 



Theorem 2.2 Let w satisfy 112. li . Then in V the mappings F, and F_ are Lipschitzfor the L2{w) 
and the ii{w) norms. Existence and uniqueness holds for l\1.2i in V H L2{w) and in V H ii{w). 

2.2 The functional CLT for converging initial data 
The time-inhomogeneous Ornstein-Uhlenbeck process 

In V, the linearization of (II. 2t around a particular solution u is the linearization of the recentered 
equation satisfied by y = g — u where g is a. generic solution for (II. 2t . It is given for t > by 

zt = 'K{ut)zt (2.2) 

where for v in V the linear operator K(u) : x K.{v)x on Cq is given by 

K{v)x{k) = aLv{k - l)^-^x{k - 1) - {aLv{kf'^ + (3)x{k) + (3x{k + 1) , k>l, (2.3) 



and is identified with its infinite matrix in the canonical basis (0, 1, 0, . . .), (0, 0, 1, ...),.. . 

A 

aLi;(l)^-^ - (aLt;(2)^-^ + /?) (3 

K(v) 



/ - {aLv{l)^-^ + P) 13 

aLv{l)^~^ - {aLv{2)^~^ +(3) (3 

aLv{2)^-^ - (aLu(3)^-i + 0) 

aLv{^)^-^ 



\ 



Let {M{k))kefi be independent real continuous centered Gaussian martingales, determined in 
law by their deterministic Doob-Meyer brackets given for t > by 

(2.4) 



{M{k))t= / {F^{u,){k) + F^{us){k)} ds . 
Jo 

The processes M = {M(k))k>o and (M) = {{M{k)))f^^-^ have values in Cq. 



Theorem 2.3 Let w satisfy ( 12. il l and uq be in V H {w). Then the Gaussian martingale M is 
square-integrable in L2{w). 
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Proof. We have E^||Mj ||^_^^^^ j = ll^^^^^) and we conclude using (I2.4t . Theorem 12.21 and 

uniform bounds in ii{w) on {us)o<s<t in function of uq given by the Gronwall Lemma. □ 



The limit equation for the fluctuations is a Gaussian perturbation of (I2.2I) . the inhomogeneous 
affine SDE given for t > by 

Zt = Zo+ [ K{us)Zs ds + Mt. (2.5) 
Jo 

A well-defined solution is called an Ornstein-Uhlenbeck process, in short OU process. We recall that 
strong (or pathwise) uniqueness implies weak uniqueness, and that li{w) C L2{w). 

Theorem 2.4 Let the sequence w satisfy ( 12.71 1. 

(a) For V in V, the operator K(w) is bounded in L2{w) with operator norm uniformly bounded in v. 

(b) Let Uo be in V H L2{w). Then in L2{w) there is a unique solution zt = e-^o ^^^"^ '^^ zq for \2.2\ 
and strong uniqueness of solutions holds for ( 12.51 1. 

(c) Let Uo be in V H ii{w). Then in L2{w) there is a unique strong solution Zt = e-fo _|_ 
JieIs^^-r)dr^M,foraandifB{\\Zo\\l^^^)) < (x^ then-E{snpt<T\\Zt\\l^i^)) < oo. 

Tightness bounds and the CLT 

The finite-horizon bounds in the following lemma will yield tightness estimates for the processes Z^ 
used in the compactness-uniqueness proof for the subsequent theorem. 

Lemma 2.5 Let w satisfy A2.ll . Let uq be in V D ii{w) and Rq be in V^. For any T > 

lim sup E ( 1 1 Zn 1 1 r / X I < oo lim sup E ( sup 1 1 Zj^ \\1 . ^ | < oo . 

We refer to Jakubowski jSl for the Skorokhod topology for non-metrizable topologies. For the 
weak topology of a reflexive Banach space, the relatively compact sets are the bounded sets for the 
norm, see Rudin fT3l Theorems 1.15 (b), 3.18, and 4.3. Hence, if i?(r) denotes the closed ball 
centered at of radius r, a set T of probability measures is tight if and only if for all e > there 
exists Te < oo such that piBir^)) > 1 — e uniformly for p in T. We state the functional CLT. 

Theorem 2.6 Let w satisfy ( 12.71 1. Consider L2{w) with its weak topology and D(M+, L2{w)) with 
the corresponding Skorokhod topology. Let uq be in V H ii{w), Rq in V^, and Z^ be given by 
(17.41 1. If {Zq)n>l converges in law to Zq and is tight, then {Z^)j\[>l converges in law to the 
unique OU process solving (12.51 1 starting at Zq and is tight. 

2.3 The functional CLT in equilibrium 

We assume the stability condition p = a/ P < 1 holds, and consider uq = u. 
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The Ornstein-Uhlenbeck process 

We set /C = K(-u) and (12.31) yields that /C : x G Cq i-^ Kx G Cg is given by 

lCx{k) =¥.{u)x{k) = l3Lp^''~\{k-l)- i^^Lp^" + I3^x{k)+ I3x{k + l), k>l, (2.6) 



identified with its infinite matrix in the canonical basis 

/ - (/3Lp^ + /3) /3 

fiLp^ - (^(5Lp^'' + /?) /3 

/3Lp^' 

V ; ; ; 






/3 



(2.7) 



Note that K, = A* where A is the generator of a sub-Markovian birth and death process. We give 
the Karlin-McGregor spectral decomposition for K, in Section IT2l to which we make a few forward 
references (it is not a resolution of the identity, see Rudin fT3l ). The potential coefficients of A are 



^ = {^{k))k>i , Ak) = L''-^piL''-L)/iL-i) ^ p-iL>^-iu{k) 



(2.8) 



and solve the detailed balance equations 7r{k + 1) = Lp^'^ T:{k) with 7r(l) = 1. The linearization of 
(ll.2t around its stable point u is the forward Kolmogorov equation for A given for t > by 



it = K.zt 



(2.9) 



which is special case of (I2.2t . Considering (I2.4t and F{u) = — F^{u) = 0, the martingale 

M = {M{k))k£N has the same law as a Cg-valued sequence B = {B{k))keN of independent centered 
Brownian motions such that B{0) =0 and for A; > 1 

v{k) := var(5i(A:)) = B{Bi{k f) = 2/3 - u{k + 1)) = 2/3p(^'-i)/(^-i) (l - /') , 

and B has diagonal infinitesimal covariance matrix diag(i)). The following result is obvious. 

Theorem 2.7 The process B is an Hilbertian Brownian motion in L2{w) if and only ifu is in ii{w). 
This is true for w = n and w = ggfor 6 > when L > 2 or for w = gefor 6 > p when L = 1. 



The Ornstein-Uhlenbeck (OU) process Z = {Z{k))kefq solves the affine SDE given for t > by 



Zt = Zo+ / KZsds + Bt 

10 



(2.10) 



which is a Brownian perturbation of (I2.9t . For L > 2, existence and uniqueness results hold under 
much weaker assumptions than (I2.lt . 
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Theorem 2.8 Let w be such that there exists c > and d > with 

< cw{k + 1) < w{k) < dp-^^^wik + 1) , k>l. 

(a) In L2{w), the operator K is bounded, the equation M.9\ has a unique solution zt = e^^ZQ 
where e'^* has a spectral representation given by ( 14. 7 I I. and there is uniqueness of solutions for the 
SDE i2.10i . The assumptions and conclusions hold for w = it and w = ggfor 9 > 0. 

(b) In addition let w be such that u is in ii{w). The SDE \2.10\ has a unique solution Zt = e^*Zo + 
Jq^^^^"^^ dBg in L2{w) further made explicit in M.2\ . This the case for w = tt and w = ge for 
9 > when L > 2 or for w = ggfor 9 > p when L = 1 . 

We use results in van Doom |l3l to prove the existence of a spectral gap, and use this fact for an 
exponential stability result inspired from Callaert and Keilson |2| Section 10. 

Theorem 2.9 (Spectral gap.) The operator fC is bounded self-adjoint in L2(vr). The least point 7 
of the spectrum of K, is such that < 7 < /3. The solution zt = e^^z^for ( 12. 91 1 in L2{n) satisfies 

\\zt\\L2{-n) < e~'^*||zo||L2(7r)- 

For L > 2 the sequence tt decays hyper-exponentially, see (I2.8I) . and (11.31) implies that the L2{it) 
norm is too strong for the CLT. Further, the mapping is not Lipschitz in V n L2 (vr) for the L2 (vr) 
norm, see Theorems 12.21 and 12.81 and their contrasting assumptions and proofs. Hence we prove 
exponential stability and (exponential) ergodicity for the OU process in appropriate spaces. 

Theorem 2.10 Let Q < 9 <1 when L>2orp<6<l when L = 1. There exists jg > and Cg < 
00 such that the solution zt = e^^zofor \2.9\ in L2{gg) satisfies \\zt\\L2{ge) — ^~'"^^^d\\zo\\L2{gg)- 

Theorem 2.11 Let w = tt or w = gg with < 9 < 1 when L > 2 or let w = gg with p < 9 < 1 
when L = 1. Any solution for the SDE 112. 1 Oi in L2{w) converges in law for large times to its unique 
invariant law (exponentially fast). This law is the law of e^^dBt which is Gaussian centered with 
covariance matrix e^*diag({i)e'^**dt made more explicit in ( 14.31 1 and \4.4\i . There is a unique 
stationary OU process solving the SDE ( 12.701 1 in L2{w). 

Global exponential stability for M.2\ . infinite-horizon and invariant law bounds, and the CLT 

We state an important global exponential stability result at u for the non-linear dynamical system. 
This is essential in the proof of the subsequent infinite-horizon bounds for the marginals of the 
processes, which yield bounds on their long time limit, the invariant law. We need uniformity over 
the state space, and results for the linearized equation M.9\ are not enough. 
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Theorem 2.12 Let p < 9 < \ and u be the solution of M.2\ starting at uq in V H L2{gg). There 
exists 70 > and Cg < oo such that \\ut — u\\i^(^gg) < e~'^^*Cg\\uQ — 'u||L2(gg). 

This does not hold in L2{tt) for L >2, else Lemma E. 13l below would also hold in L2(vr), which 
would contradict (ll.3t . Theorem 3.6 in Mitzenmacher 1121 states a related result for some weighted 
ii norms obtained by potential function techniques. 

Lemma 2.13 Let p < 6 < 1 when L>2orp<6<l when L = 1. Then 

limsupE ( < ^ limsupsupE ( , J < oo 

and under the invariant laws limsup^_^oo (H ■^0^11^2(99)) ^ 

Our main result is the functional CLT in equilibrium, obtained with a compactness-uniqueness 
method using tightness of the invariant laws (based on Lemma l2.13t and Theorems l2.6l and l2.lll 

Theorem 2.14 Let the networks of size N be in equilibrium. For L > 2 consider L2{gp) with 
its weak topology and ]D)(M+, L2{gp)) with the corresponding Skorokhod topology. Then {Z'^)n>l 
converges in law to the unique stationary OU process solving the SDE ( 12. 7 OH . in particular {Zq)j\^> l 
converges in law to the invariant law for this process (see Theorem \2.11l . For L = 1 the same result 
holds in L2{ge)for p < < 1. 

3 The proofs for converging initial conditions 

3.1 Existence and uniqueness results 

Proof of Theorem 12.21 (refined existence result for il.2i ) 

We give the proof for L2{w), the proof for li{w) being similar. The assumption (I2.lt and the identity 

-y^ = {x- y){x^'^ + x^~'^y H h y^'^) yield 

{u{k - l)^ - v{k - l)^fw{k)-^ < {u{k - 1) - v{k - l)f L^dw{k - 1)"^ , 
{u{k)^ - v{k)^)^ w{k)-^ < {u{k) - v{k)f L^w{k)-\ 
{u{k + 1) - v{k + l)f w{k)-^ < {u{k + 1) - v{k + l)f c~^w{k + 1)-^ , 

hence we have the Lipschitz bounds — < 2a^L'^{d + l)||n — 'y||^2(^) and 

||-F_(m) — < 2/3^(c~-'- + l)||u — Existence and uniqueness follows by a 

classical Cauchy-Lipschitz method. 
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The derivation of the Ornstein-Uhlenbeck process 

Let {x)k = x{x — 1) • • • (x — /c + 1) for x G M (the falling factorial of degree k G N). Considering 
(ll.lt . let the mappings and with values in Cg be given for v in cq by 

Ffwm = J^""^-'>>V"^"""' - '"'(') = F!;iv)-F.(v)^ ,3.1) 

The process is Markov on V^, and when in state r has jumps in its k-th coordinate, A; > 1, of 
size 1/N at rate NF^{r){k) and size -1/iV at rate NF_{r){k). 



Lemma 3.1 Let Rq be in V^, u solve ( li.2l) starting at uq in V, and be given by ( 17.41) . Then 

Z[^ = Z^+ f' y/N{F^{Rf)-F{us))ds + Mt^ (3.2) 
Jo 

defines an independent family of square-integrable martingales = {M^ {k))k^f^ independent 
of Zq with Doob-Meyer brackets given by 

(M^(fc))^= f {F^{Rf){k)+F_{Rf){k)}ds. (3.3) 
Jo 

Proof. This follows from a classical application of the Dynkin formula. □ 

The first lemma below shows that it is indifferent to choose the L queues with or without replace- 
ment at this level of precision. The second one is a linearization formula. 

Lemma 3.2 For N > L > 1 and a in M we have 



A\a):='-^-a^ = Y.{a-iya^-= E 



and A^{a) = N-^0{a), uniformly for < a < 1, and A'^{k/N) < 0/or k = 0, 1, ... , N. 

Proof. We develop = U^'q = \{^~^ {a + {a- to obtain the identity for 

A^{a) which is clearly N~^0{a), uniformly for < a < 1. For a = k/N, Hf'^ is com- 

posed of terms bounded by a or contains a term equal to and cannot exceed a^. □ 

Lemma 3.3 For L > 1 and a and h in R we have 



B{a, h) := (a + h)^ - - La^^^h = ^ (^^ 

,_9 V * / 



1=2 

with B{a, h) = Ofor L = I and B{a, h) = for L = 2. For L > 2 we have < B{a, h) < 
+ {2^ - L - 2) ah'^for a and a + hin [0, 1]. 
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Proof. The identity is Newton's binomial formula. A convexity argument yields B{a, h) > 0. For a 
mda + h in [0, 1], B{a, h) < + Ylf=2 = + {2^ - L - 2) ah^ . □ 

Let V be in V and x in Cg. Considering (ll.lt . (13.11) and Lemma lT2l let : V ^ Cg be given by 

{v){k) = aA^ {v{k - I)) - aA^ {v{k)) , k>l, (3.4) 
and considering dl.ll) . (12.31) and Lemma l33] let i/ : V x cj] ^ Cg be given by 

H{v,x){k)=aB{v{k-l),x{k-l))-aB{v{k),x{k)), k>l, (3.5) 
so that for f + X in V 

F^ = F + G^, F{v + x)-Fiv)=K{v)x + H{v,x), (3.6) 
and we derive the limit equation (12.51) and (I2.4t for the fluctuations from (13. 2t and (13. 31) . 

Proof of Theorem 12.41 (existence and uniqueness for the OU process) 

Considering (12. 3t . < 1, convexity bounds, and (I2.lt . we have 

W^i^Mhiw) < 2(aL + /3) ^ (aLx(A: - lfdw{k - 1)"^ + {aL + /3)x(A;)2u;(A;)-^ 

fc>i 

+ /3x(/£ + l)2c"^w(A: + l)-^) 
< 2{aL + (3){aL{d + 1) + (3{c-' + l))||x||i^(^) 

and (a) and (b) follow. For Uq in V Pi ii{w) the martingale M is square-integrable in L2{w). If 
E^IIZqIII^^^^^ < oo then the formula for Z is well-defined, solves the SDE, and the Gronwall 
Lemma yields E^supj<y ||^_^^^^^ < oo. Else for any e > there is < oo such that 
P(||Zg||/^2(«>) and a localization procedure using pathwise uniqueness yields existence. 

3.2 The proof of the CLT 

Proof for Lemma 1131 (finite-horizon bounds) 

Using and 

= + + y/N f G^{R^) ds + f y/N {F{R^) - F{u,)) ds (3.7) 
where LemmalTlvields G^{R^){k) = N-^0{R^{k - 1) + R^ {k)) and considering dO 

l|G"^(^f)IL.H=^"^0(ll^flL.(.))- (3-8) 

We have 
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and since F+, F_ and F are Lipschitz (Theorem 12.21) the Gronwall Lemma yields that for some 
Kt < CO we have - \\uo\\l2{w) ^nd 

sup llzf |L , , <Kt( \\z^L , , + N'^/'^KtWuqWj ,s + sup ||m/^|L , J. 

We conclude using the Doob inequahty, (I3.3t . (I3.6t . (I3.8t . (I3.9I) . and 

||i^+(i^f ) + ^~«)IL,(^) < ll^f IL,(.) • (3-10) 

Tightness for the process 

Lemma 3.4 Let w satisfy 112. li . and consider L2{w) with its weak topology and D(M-|-, L2{w)) with 
the corresponding Skorokhod topology. Let uq be in V H £i{w) and Rq in V^, and be given by 
AL4\l . If (Z^) 7v> L is tight then (Z^) n>l is tight and its limit points are continuous. 

Proof. For e > let < cxo be such that ^{Z^ ^ B{rs)) > 1 — e for N > 1 (see the discussion 
prior to Theorem l2.6t . Let R^'^ be equal to R^ on {Z^ G B{ri;)} and such that Z^'^ is uniformly 
bounded in L2{w) on {Z^ B{r^)}. Then Z^'^ is uniformly bounded in L2{w) and we may 
use a coupling argument to construct Z^'^ and Z^ coinciding on {Z^ G B{r^)}. Hence to prove 
tightness of (Z^)7v>l we may restrict our attention to {Zq)]\j>l uniformly bounded in L2{w), for 
which we may use Lemma l231 

The compact subsets of L2{w) are Polish, a fact yielding tightness criteria. We deduce from 
Theorems 4.6 and 3.1 in Jakubowski |j8|, which considers completely regular Hausdorff spaces (Ty- 
chonoff spaces) of which L2{w) with its weak topology is an example, that {Z^) Ar> ^ is tight if 

1. For each T > and e > there is a bounded subset Kr^e of L2{w) such that for > L we 
have P(Z^ G D([0, T], Kt,^)) >l-e. 

2. For each d> 1, the d-dimensional processes (Z^(l), . . . , Z^ {d))N>L are tight. 

Lemma 1231 and the Markov inequality yield condition 1. We use (13. 7t (see (13. 2t and (13. 6t ). 
and (13. 3t and (13. 6t . and the bounds (13. 8t . M.9\ and (I3.10t . The bounds in Lemma 1231 and the fact 
that Z^{k) has jumps of size l/\/iV = o{N) classically imply that the above finite-dimensional 
processes are tight and have continuous limit points, see for instance Ethier-Kurtz |4|| Theorem 4.1 
p. 354 or Joffe-Metivier (9) Proposition 3.2.3 and their proofs. □ 

Proof of Theorem nil (the functional CLT) 

Lemma l3!4l implies that from any subsequence of Z^ we may extract a further subsequence which 
converges to some with continuous sample paths. Necessarily Zq° has same law as Zq. In (I3.7t 
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we have considering (I3.6t that 

VN{F{R^){k) - F{us){k)) = K(n,)Zf + VnH{us, N-^'^Z^) . (3.11) 

We use the bounds (I3.8I) . i3.9i and (I3.10t . the uniform bounds in Lemma l231 and additionally (I3.5t 
and Lemma l331 We deduce by a martingale characterization that Z°° has the law of the OU process 
unique solution for (I2.5t in L2{w) starting at Z^, see Theorem I2.4t the drift vector is given by the 



limit for (I3.2I) and (I3.7I) considering (I3.1 II) . and the martingale bracket by the limit for (13.31) . See for 
instance Ethier-Kurtz Theorem 4.1 p. 354 or Joffe-Metivier j9|| Theorem 3.3.1 and their proofs 
for details. Thus, this law is the unique accumulation point for the relatively compact sequence of 
laws of {Z^)]\j>i, which must then converge to it, proving Theorem 12. 61 

4 The properties of AC = K({t) 

4. 1 Proof of Theorem 12.81 (existence and uniqueness results) 

Considering (12. 6t and convexity bounds we have 

W^^U^) = E ( - 1) - (V + iHk) + zik + 1)) ' w{k)-' 

k>l 



k>l k>l 

+ Y,4kfw{ky^+Y,z{f^ + ^?'w{f^)~^) 

k>l k>l ^ 

< f(2L + 2) ( Ld^z{k - \fw{k - 1)"^ + ( + 1) A^f^ik)-^ 

k>2 k>l 



+ c-^^z{k + lfw{k + l)-A 

k>l ^ 



< f3\2L + 2) ( ^ ^ ^-1 ^ p||2^^^^ _ 
The Gronwall Lemma yields uniqueness. For > 1 we have 

(Lp^)"' TTik + 1) < nik) = [Lp'^'y^ 7r{k + 1) < (L-^/p"'''') vr(fc + 1) , 

When i? is an Hilbertian Brownian motion, the formula for Z yields a well-defined solution. 
4.2 A related birth and death process, and the spectral decomposition 



Considering (I2.7t . A = IC* is the infinitesimal generator of the sub-Markovian birth and death 
process on the irreducible class (1, 2, . . .) with birth rates = PLp and death rates p.^ = (3 for 
k>l (killed at rate pi = (3 at state 1). The process is well-defined since the rates are bounded. 
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Karlin and McGregor (TOIITTJ give a spectral decomposition for such processes, used by Callaert 
and Keilson EE) and van Doorn |3| to study exponential ergodicity properties. The state space in 
these works is (0, 1, 2, . . .), possibly extended by an absorbing barrier or graveyard state at —1. We 
consider (1, 2, . . .) and adapt their notations to this simple shift. 

The potential coefficients (| 10 1 eq. (2.2), [3] eq. (2.10)) are given by 

and solve the detailed balance equations /ifc+i7r(A; + 1) = Xk7r(k) with 7r(l) = 1, see ( I2.8t . 

The equation AQ{x) = —xQ{x) for an eigenvector Q{x) = {Qn{x))n>i of eigenvalue —x 
yields \iQ2{x) = (Ai + /Ui - x)Qi{x) and \nQn+i{x) = (A„ + fin - x)Qn{x) - HnQn-i{x) for 
n > 2. With the natural convention Qo = and normalizing choice Qi = 1, we obtain inductively 
Qn as the polynomial of degree n — 1 satisfying the recurrence relation 

-xQn{x) = f3Qn-i{x) - {pLp^" + f3) Qn{x) + /3L/"g„+i(x) , n > 1 , 

corresponding to 1 10 1 eq. (2. 1) and |3|| eq. (2. 15). Such a sequence of polynomials is orthogonal with 
respect to a probability measure -0 on and, for i,j>l with i ^ j, Qi{x)'^ ijj{dx) = vr(i)^^ 
and /g°° Qi{x)Qj{x) 'ip{dx) = or in matrix notation Q{x)Q{x)* i/j{dx) = diag(7r^-'^). 

Let Pt = {pt{hj))i,j>i denote the sub-stochastic transition matrix for A. The adjoint matrix 
Pf is the fundamental solution for the forward equation zt = A*zt = ICzt given in (12.91) . The 
representation formula of Karlin and McGregor fTOlfTTl . see (1.2) and (2.18) in [3 1, yields 

/•OD 

e'C* = P; = {p*tihj)kj>i , Plihj) =PtU,i) = <i) / e-^'Q,{x)Q,{x)^p{dx) , (4.1) 

Jo 

or in matrix notation e^* = diag(7r) e~^*'Q{x)Q{x)* t(j{dx). 

The probability measure tp is called the spectral measure, its support S is called the spectrum, 
and we set 7 = minS". The OU process in Theorem l2. 81 (b) and its invariant law and its covariance 
matrix in Theorems l2.1 ll and l2.14l can be written 

Zt = diag(7r) J e-^'*Q(x)* (^Zq + e^^ dBs^ Q{x) i>{dx) , (4.2) 
e^' dBt = dmgi'K) j \Q{xY j e-""' dBA Q{x) ij{dx) , (4.3) 
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e^*diag({))e'=** dt = diag(7r) [ Q(^)*diag(^)Q(j/) g^^^g^^)* ^{dx)i>{dy) diag(^). (4.4) 
Js2 x + y 

4.3 The spectral gap, exponential stability, and ergodicity 

Proof of Theorem 12.91 (spectral gap and exponential stability in the self-adjoint case) 



The potential coefficients ('/r(A;))fc>i solve the detailed balance equations for A and hence JC = A* is 
self-adjoint in L2(vr). For the spectral gap, we follow Van Doorn [3 1, Section 2.3. The orthogonality 
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properties imply that Qn has n — 1 distinct zeros < Xni < ■ ■ ■ < Xn,n-i such that Xn+i,i < Xn,i < 
Xn+i,i+i for 1 < i < n - 1. Hence = Yin\n^ooXn,i > exists, < fj+i, and a = limj^oo 
exists in [0, oo]. Theorem 5.1 in |3l estabhshes that 7 > if and only if a > and Theorem 5.3 (i) 
in that a = {yj\\mk Xk — \/limfc Hk) ^ = /3 > 0. (Theorem 3.3 in states that 7 = ^1 < cr, but 
estimating ^1 is impractical.) 

For the exponential stability, we have ||^t||^2{7r) ~ (^'^*'^0) e^* -20)2,2(77) ^® ^^^^ 
self-adjoint in L2(7r) and the spectral representation (14.11) yield 



(e'^*zo,e'^*zo)^^(^) = {zo,e''^' zq) ^^^^^ = J ^e-^^'z*oQ{x)Q{xr zoiP{dx) 

< / z*oQix)Qixrzo^{dx) = e'^^* (^o,^o)l2W ' 



s 

Proof of Theorem IZlOl (exponential stability, non self-adjoint case) 

It is similar to and simpler than the proof for Theorem 12. 121 to which Section |5l is devoted, and we 
postpone the proof until the end of that section. 

Proof of Theorem |2.1 II (ergodicitv for the OU process) 

We use the uniqueness result and explicit formula in Theorem l2.8l and Theorem l2.9l or l2!T0l 

5 Exponential stability for the nonlinear system 
5.1 Some comparison results 

Considering (13. 6t with /C = K(tt) and F{u) = 0, if n solves (ll.2t in V then y = u — u solves the 
recentered equation given by yt{k) = F{u + y) = K,yt{k) + H{u, yt){k) or 

yt{k) = PLp^^'^ytik - 1) + a5(u(fe - l),yt{k - 1)) 

- (^pLp^'ytik)+aBiu{k),ytik))+pytik)) +Pytik + l), k>l. (5.1) 

If no is in V n £1 then n is in V H ^1 and hence y is in £^ and for A; > 1 we have 

Uk) + Uk + 1) + • • • = PLp^'^'ytik - 1) + aB{uik - l),yt{k - 1)) - Pyt{k) . (5.2) 

If y solves (I5.lt starting at y^ such that yo + is in V, then u = y + u solves (ll.2t in V starting at 
^0 = 2/0 + u- Then —u<y< 1 — u and — 1 < y < 1. For yo + uinV H iiwe have y in i^. 



Lemma 5.1 Let u and v be two solutions for t\L2} in V such that uq < vq. Then ut < vtfor t > 0. 
Let yo + ube in V and y solve (15.71 1. Ifyo > then yt > and if y^ < then yt < Ofor t > 0. 
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Proof. Lemma 6 in 1151 yields the result for (ll.2t (the proof written for L = 2 is valid for L > 1). 
The result for (I5.lt follows by considering u = y + u and u which solve (II. 2t . □ 

We compare solutions of the nonlinear equation (I5.lt and of certain linear equations. 

Lemma 5.2 Let A be the generator of the sub-Markovian birth and death process with birth rate 
Afc > and death rate P at k > 1. Let sup^ < oo. The linear operator x ^ A.*x given by 

A*x{k) = Xk^ix{k - 1) - (Afc + P)x{k) + (3x{k + 1) , k>l, 

is bounded in ^5- There exists a unique z = {zt}t>o given by zt = e"^ *zo solving the forward 
Kolmogorov equation z = A* z in It is such that if zq > then zt > and if zq < then zt < 0, 

and zt{k) + zt{k + 1) H = Xk-iZt{k - 1) - Pzt{k) for k > 1. 

Proof. The operator norm in £^ of A* is bounded by 2(supfc Afc+/?), hence existence and uniqueness. 
Uniqueness and linearity imply that if zq = then zt = and else if zq > then Zf||zo|lr^ "^^e 
instantaneous law of the process starting at zolkoIlT^ hence zt > 0. If < then —z solves the 
equation starting at —zq > and hence —zt>0. The last result is obtained by summation. □ 

Lemma 5.3 Let L > 2 and y = {yt)t>o solve ( 15. 7 1 with yQ + u inV D i\. Under the assumptions of 
Lemma U^ let z = {zt)t>o solve z = A*z in Let h = {ht)t>o be given in £^ by 

h{k) = z{k) +z{k + l) + -- - - {y{k) +y{k + l) +■■■), k>l. 

(a) Let Afe > (5Lp^'' + a(l + (2^ - L - 2) u{k)) for k>l,yQ>Q, and ho > 0. Then ht > Ofor 
t > 0. 

(b) Let Afc > PLp^'^for k > 1, yo < 0, and Hq < 0. Then ht < Ofor t > 0. 

Proof. We prove (a). For e > let A* correspond to A|, = Afc + £■ The operator norm in 
of A* — A* is bounded by 2e, hence lime_*o e-^^^zo = zt in and we may assume that A^ > 
f3Lp^'' + a(l + (2^ — L — 2) u{k)^ for k > 1. Since zt = e-^'^zo depends continuously on zq in 
we may assume /iq > 0. Let r = mi{t > : {A; > 1 : = 0} 7^ 0} be the first time when 

h{k) = for some A; > 1. We have r > 0. 

The result (a) holds if r = 00. If r / 00, Lemma ls!2l and (15. 2t yield 

hrik) = Xk-Wrik - 1) - (3Lp^'"\r{k - 1) - aB{u{k - - 1)) 

+ Xk-i{zr{k - 1) - yr{k - 1)) - P{zr{k) - yr{k)) . 
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Lemma ISlTI yields y > and Lemma 1331 and y < I yield 

B{u{k - 1), y{k - 1)) < y{k - 1)^ + (2^ - L - 2) u(A; - l)y{k - if 
< (1 + (2^ - L - 2) u{k - 1)) y(A; - 1) , 

^ r fc — 1 

hence Xk-iy{k — 1) — (3Lp y{k — 1) — ciiB{u{k — l),y{k — 1)) > with equality only when 
y{k - 1) = 0. For /c in Z = {fc > 1 : hr{k) = 0} / we have 

Zr{k - 1) - yr{k - 1) = hr{k - 1) > , Zr{k) - yr{k) = -hr{k + 1) < , 

hence hrik) > with equality if only if /c — 1 is in Z U {0} and /c + 1 is in Moreover ht{k) > 
for t < r and hr{k) = imply hrik) < 0. Hence hT-{k) = 0, and the above signs and equality 
cases yield that Zr{k — 1) = yr(k — 1) = and /c — 1 is in Z U {0} and A; + 1 is in By induction 
ZT-{i) = yrii) = for i > 1 which implies zt = yt = for t > r, and the proof of (a) is complete. 

The proof for (b) is similar and involves obvious changes of sign. The assumption > j3Lp^^ 
suffices to conclude since B{u{k — 1), y{k — 1)) > (Lemma l3.3t and the non-linearity "pushes" 
in the right direction. □ 

Lemma 5.4 For any < < 1 there exists Kg < oo such that for x in L2{ge) C ^5 

\\{x{k) +x{k + !) + ■■ •)fc>i 11^,(5,) < ^elkllL2(5fl) • 

Proof. Using a classical convexity inequality 

^(x(A;) + x(fc + l) + ---)'^~'^ 
k>l 

< ^n{x{kf + x{k + lf ^ h x(A; + n - 2)^ + {x{k + n - 1) + x{k + n) ^ f) 6'^ 

k>l 

< n(i + + • • • + x{kfe-^ + n r-^ Y{x{k) + x{k + i) + --- fe-^ 

k>l k>l 

and we take n large enough that n^'^"^ < I and = n{l + ^ h 6'"'"2) (1 - nfl""^)"^ □ 

5.2 Proofs of the exponential stability results 
Proof of Theorem ITTH for L > 2 

If no is in V n L2{ge), then so are = m.m{uo,u} and Uq = maxjuo, tt} and hence the corre- 
sponding solutions and ti+ for (ll.2t . see Theorem l2.2l Lemma lSlTl vields that < ut < uf and 

< u < uf for t > 0. Then 

y = u — u , = — u > , = — u < , 
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solve the recentered equation ( I5.lt . and termwise 

|yo| = max{y(J", -y^} , \yt\ < max{y^ , -y^} , t > . (5.3) 
We consider the birth and death process with generator A defined in Lemma ls!2l with 

Afc = max|/?L/j^' +a(l + {2^ - L - 2) u{k)) ,/3(9| , k>l, 

which satisfies the assumptions of Lemma 1531 (a) and (b). We reproduce the spectral study in Sec- 
tion l4.2l and the proof of Theorem l2.9l in Section l43] f or A, corresponding objects being denoted with 
a hat. For p < 6* < 1 we have a < (59 and hence is equivalent to (39 for large k, Theorems 5.1 
and 5.3 (i) in Q yield that < 7 < a = {y/j39 - y/^f = (3 {l - ^f, and if z solves i = A*z 

then ll^t ||l2{*) ^ 6~^*||2:o||l2(*) for ^ ^ 0- Moreover 

k-l 

9^-^ < TT{k) = 9^-^ Yl max {^"^Lp^' + 9-^p{l + (2^ - L - 2) u{k)) , l} 

i=l 

and the product converges, hence 7r{k) = 0{9^) and 9^^ = 0{'K{k)) and Lemma lTTl vields that there 
exists c > and d > such that c^^\\ ■ H^jW — II " Ili2(ge) — '^ll " ll-L2{^-)- Hence for t > 

||^t||L2(ge) ^ t^ll^tlli2(*) ^ e~^*t?lko||L2{,}) < 6-^*0^112011^2(3,) . 

Hence if z+ solves z+ = A*z^ starting at Zq = y^ > then Lemmas 15.31 (a) and l5.4l vield 

\\yt\\L2{ge) < \\iytik)+ytik + 'i-) + ---)k>l\\L2(gg) 

< \\i4ik)+4ik + l) + ---)k>l\\L,(ge) 

< KeWzth^ige) < e^^*c(iKe||y+||i2(gg) 

and similarly if z~ solves z~ = A*z^ starting at z^ = y^ < then Lemmas 15.31 (b) and l5.4l vield 
||yt~llL2{39) ^ e-T'*cdK0||yo We set 7^ = 7 and Ce = cdKe. Considering (|53J, 

II l|2 ^ II +||2 I II -||2 ^ -2'-fgt/^2 ( II +||2 I II -||2 \ 

\\yt\\L2{ge) - \\yt \\L2{ge) + H^^t \\L2{ge) - ^ yWo \\L2{go) + H^O \\L2{ge)) 

and we complete the proof by remarking that for A; > 1, either ?/o'(fe) = yo(^) and y^ik) = or 

Voik) = yo{k) andy+{k) = 0, and hence \\y^\\l^^g^) + lly(7llL{5e) = llyoHiaCse)- 

Proof of Theorem ITTol and of Theorem ITTH for L = 1 

The linearization ( I2.9t of Equation (II. 2t is obtained by replacing B and in Equation (I5.lt by 
and coincides with Equation (I5.lt for L = 1. Likewise, the equation for ( I2.9t corresponding to (15. 2t 
is obtained by omitting the terms aB{u{k — l),yt{k — 1)). We obtain a result for Equation ( I2.9t 
corresponding to Lemma 1531 (a) and (b) under the sole assumption > jSLp^'' for k > 1. The 
proof proceeds as for Theorem l2.12l for L > 2 with the difference that A^ = max {(3Lp^\(39}. We 
have Afc equal to (39 for large A; for < < 1 when L > 2 and for p < ^ < 1 when L = 1. 
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6 Tightness estimates and the functional CLT in equilibrium 

6.1 Proof of Lemma l2lT3l (infinite horizon and invariant law bounds) 

Let Uh{v) be the solution of (ll.2t at time h > with initial value v in V. For to ^ let = 

N {R^^^^ - [//.«)) . Then Z,^^^ = Z.^^ + VN ([/,«) - u) and Theorem m yields 

The conditional law of {Z^ fJh>o given = r is the law of Z^ started with Rq = uq = r, in 

■N 

to,' 



particular with Zq = Z^ o ~ 0. We reason as in (I3.7l) - (l3.10t except that the bound (I3.9I) becomes 



WnN II ^ IL-^II _L Ar-1/2 II || 

and we use (I6.1I) and obtain that for some Kt < oo 

which combined with (I6.1I) yields that for some Lt < oo we have for < /i < T 

E < + 2(^tA^-^ + e--«^)^C,^ E • (6.2) 

We fix T large enough for 8e-2Te'^C| < e < 1. Uniformly for N > Kto^oT ^ for m G N 

yN ||2 \ ^ T , /ll vAf l|2 



and by induction 

7N ||2 \ / , ^rm:, /"ll^ATIlS \ ^ , ^/'||^Ar||2 



and (I6.2t also yields 
hence the infinite horizon bound 

Ergodicity and the Fatou Lemma yield that for Z^ distributed according to the invariant law 



<lT>irfE(^||Zn|,,(,J <supE(^||Z, 



and the invariant law bound follows if we show that we can choose R^ in such that 

7Af ||2 
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limsupE(||Zo^||'^,^^_J <oo. (6.3) 



For this we consider L > 2, the case L = I being similar, and Rq given for k >Ohy R^{k) = 
iN-^ with l<i<N such that < u{k) - iN'^ < 2-^N-\ For x > and < ?/ < 1 

y = p{L^'l)/(L-l) 4^ X = log (1 + (L - 1) log y/ log p) I log L 

4^ ^-^ = (1 + (L - l)logy/logp)''°s^/^°s^ 

hence for z{N) = infjfc > 1 : R^(k) = O} we have z{N) = mf{k > 1 : u{k) < 2-'^N-'^} = 
mf{ken:k> log (1 + (L - 1) log {2-^N-'^) /logp) /logL}. Then 

z{N)-l 

k=l k>z{N) 

with 

N J2 {Roik)-u{k)) <2-^N~' ^ =0(jV"i(logiV)-^°g^/'°s^) 
fc=i 

and for large enough (and hence z(N)) 

k>z{N) j>0 

hence (I6.3I) holds and the proof is complete. 

6.2 The functional CLT: Proof of TheoremEU 

Lemma 12.131 and the Markov inequality imply that in equilibrium {Zq)n>l is tight for the weak 
topology of L2{gp), for which all bounded sets are relatively compact. Consider a subsequence. We 
can extract a further subsequence along which {Zq)n> l converges in law to some square-integrable 
in L2{gp), and Theorem 12.61 vields that along the further subsequence {Z^)n>l converges in 
law to the OU process unique solution for (I2.10t in L2{gp) starting at Z^. 

The limit in law of a sequence of stationary processes is stationary (Ethier-Kurtz f4l p. 131, 
Lemma 7.7 and Theorem 7.8). Hence the law of Z°° is determined as the unique law of the stationary 
OU process given by (I2.10I) . see Theorem 12.111 From every subsequence we can extract a further 
subsequence converging in law to Z°°, hence liniAr^oo = Z°° in law. 
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